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Abstract 

We study the counting statistics of charge transport in the Anderson impurity model (AIM) 
employing both Keldysh perturbation theory in a Fermi liquid picture and the Bethe ansatz. In 
the Fermi liquid approach, the object of our principal interest is the generating function for the 
cumulants of the charge current distribution. We derive an exact analytic formula relating the 
full counting statistic (FCS) generating function to the self-energy of the system in the presence 
of a measuring field. We first check that our approach reproduces correctly known results in 
simple limits, like the FCS of the resonant level system (AIM without Coulomb interaction). 
We then proceed to study the FCS for the AIM perturbatively in the Coulomb interaction. By 
comparing this perturbative analysis with a strong coupling expansion, we arrive at a conjecture 
for an expression for the FCS generating function at 0{V^) (V is the voltage across the impurity) 
valid at all orders in the interaction. 

In the second part of the article, we examine a Bethe ansatz analysis of the current noise for the 
AIM. Unlike the Fermi liquid approach, here the goal is to obtain qualitative, not quantitative, 
results for a wider range of voltages both in and out of a magnetic field. Particularly notable are 
finite field results showing a double peaked structure in the current noise for voltages satisfying 
eV ~ fJ-sH. This double peaked structure is the "smoking gun" of Kondo physics in the current 
noise and is directly analogous to the single peak structure predicted for the differential conductance 
of the AIM. 

PACS numbers: 
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I. INTRODUCTION 

The subject of counting statistics is rooted in the historical paper by Schottky [l| where 
the measurements of charge noise have been carried out and interpreted as the basis for 
determining the elementary charge e of the current carriers: electrons. Contemporary trans- 
port experiments are being performed on nano-structures, usually involving two electron 
reservoirs (left and right) and a central constriction j^. The mean electric current, or linear 
conductance is well understood in terms of scattering theory and, for a single conducting 
channel, is given by the Landauer formula: 

Go = , (1) 

where Tq is the transmission coefficient and factor 2 stems from the electron spin. However, 
due to the quantum nature of the problem, the current is bound to fluctuate. In particular, 
this gives rise to interesting noise (the second moment of the current distribution) properties 
extensively discussed in the literature \4\. 

With the third moment of the current distribution now available experimentally [5|, it is 
natural to widen the question to the full current distribution function or the full counting 
statistics (FCS). One way to formulate this question is to ask what is the probability P{Q) 
that charge Q will be transmitted through the system during the waiting time T and for a 
given bias voltage V. As electrons are discrete particles, a naive guess at P{Q) would be 
the Poisson's distribution: 

PiQ) = , (2) 

where (Q) = GqVT. For simplicity, we set e = = 1 in what follows. The electrons, 
however, are not only discrete particles but also quantum particles obeying Fermi-Dirac 
statistics. Due to the Pauli principle the electrons will tunnel 'one by one'. So, given the 
'number of attempts', = VT/tt, one would expect the total probability be proportional 
to the probability of successes Tq as well as the probability of failures (1 — Tq)^~^ . The 
resulting probability distribution is binomial: 



P{Q) 



Q 

N 



To«(l-To)^-«, (3) 



where the binomial coefficient in front simply follows from the normalization: PiQ) = 1- 
Note that the binomial distribution is a clear signature of Fermi statistics; indeed, the re- 



spective probability distribution for bosons is the inverse binomial 6j. In practice, it is more 
convenient to work with the generating function xW = Xlg ^^'^''^-^(Q)' where the Fourier 
transform variable A is called the 'counting field' (see below). The irreducible moments of 

Q 

the charge distribution immediately follow [7|: 

lnx(A) = f;((Q"))^. (4) 

n=l 

The generating function for the binomial distribution is simply: 

Xbinomial=[l + ro(e'^-l)]^. (5) 

From this equation one easily recovers the Landauer formula, (Q) = NTq, the well known 
expression for the shot noise {{Q"^)) = NTq{1 — Tq), and obtains the following expression for 
the third moment: {{Q^)) = iVTo(l — To)(l — 2To). Note that for low transmission (Tq 0), 
the statistics reverts to Poissonian, while for perfect transmission (To 1), there are no 
current fluctuations and xW = «AA^. The physics described so far has been understood in 
the seminal paper by Levitov and Lesovik 8] (see also j^), where they derive a more general 
formula for the generating function 

lnxo(A;y;{T(c.)}) = 2r J ^ln{l+T(c.) (6) 

X [uLil - n^)(e*^ - 1) + nn{l - nL){e-'^ - 1)] } , 

which is valid for finite voltage, temperatures, and allows for the energy dependent trans- 
mission coefficient. Here nL/R{uj) = npiuj =F V/2) are the thermal electron distributions in 
the left and right leads and np{uj) is the Fermi function. Clearly Tq is T{uj) at the Fermi 
energy, set at u = 0. Schonhammer has recently re-examined this formula, Eq.([n]), by an 



io|. 



alternative method, and found that it is correct 

The discussion so far has focused upon non-interacting electrons. But while counting 
statistics for non-interacting electrons is by now comprehensively understood, the same 
cannot be said when the electrons interact with each other and with the substrate. Conse- 
quently the understanding of the interaction effects on the FCS has become an important 
issue. There has been many papers on the subject in recent years with many interesting 
yet miscellaneous results. Certainly no general paradigm as to how interactions should af- 
fect the statistics has as yet emerged. Misunderstandings dominating the subject only a 



4 



short time ago are well illustrated by the following example. The generating function xW 



for Matveev's Coulomb blockade setup 



ll| (equivalent to the g = 1/2 Kane and Fisher 



problem 12L which is in turn equivalent to the a = 1/2 dissipation problem first so 



by Guinea 
and 



131]) has been calculated independently by three different methods in 14j . 



ved 



id, 



161 ] with seemingly very different results. It was only understood later that all three 



results are indeed correct and represent one and the same function (see Appendix C to Ref. 



171]). Moreover the distribution in question turned out both to be simple and to represent 



a particular case of Eqn. with a specific choice of the transmission coefficient, T{u). 
This lack of a coherent picture of the FCS in strongly correlated systems is, we believe, 
simply explained. As is illustrated by the prominence of the Fermi and Luttinger liquid 
paradigms, it is accepted that in the condensed matter it is the low-energy physics which is 
universal. The FCS is no exception. On an energy scale set by the bare parameters, it is 
therefore the low-temperature, low-voltage expansions of xW where universal results are to 
be found. The high-voltage (temperature) distributions may be enormously fascinating but 
are destined to remain model dependent. 

In this paper we collect together a number of such universal results, presented in two parts. 
In the first part. Sections II and III, we study the generating function, x, using Keldysh 
perturbation theory in a Fermi liquid approach. In Section II, we introduce the Keldysh 
method for calculating the statistics. In the process we establish an exact relationship 
between the generating function and the self energy. In Section III, we study the FCS for 
the AIM both in perturbation theory and in the strong coupling limit. By comparing the 
two we propose a conjecture for x at low voltages, i.e. 0{V^), but valid at all orders of the 
interaction. 

In the second part of the paper. Sections IV through VI, we switch tacts and instead 
employ a Bethe ansatz analysis of the current moments for the AIM. We however limit 
ourselves to exploring the behavior of the current and the current noise in the AIM's Kondo 
regime. Our results for these quantities differ from those of the first section of this paper. 
Here the focus is on their qualitative not quantitative features but over a larger range of 
voltages (though still much smaller than any bare energy scale) and for finite magnetic 
fields. In Section IV we review the Bethe Ansatz method for calculating the current and 
noise. We then present results for zero magnetic field in Section V where, in addition, a 
comparison is made to Fermi liquid calculations. In the final section. Section VI, we consider 
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the properties of the noise in finite magnetic fields. Here is found the most significant result 
of the second part of the paper. We are argue that in the vicinity of voltages commensurate 
with the magnetic field, the current noise should see a double humped enhancement. This 
enhancement is the analog of that seen in the current when eV ~ fJ^sH and so represents a 
'smoking gun' 46|| of Kondo physics. 



II. KELDYSH METHOD FOR THE CALCULATION OF CURRENT STATIS- 
TICS: GENERAL RESULTS 



The calculation of the charge statistics is usually accomplished by coupling the system 
to a 'measuring device'. In the original gedanken experiment by Levitov and Lesovik it 
is a fictitious spin-1/2 galvanometer coupled to the current 9(]. The transmitted char ge i s 
then proportional to the change of the spin phase. As has been shown by Nazarov 18| . 
the counting of charge can in general be done by coupling the system to a fictitious field 
and calculating the non-linear response, so leading to the same results. In fact the standard 
quantum mechanical formula, P{Q) = {S{Q — Q)), can also be used provided that the central 
region is initially decoupled from the leads. 



According to 



19| the generating function is given by the following average. 



X(A) = ( Tc exp 



-z / Tx{t)dt 



(7) 



where C is the Keldysh contour [20|], A(t) is the measuring field which is non-zero only during 
the measuring time T: A(t) = X6{t)6{T — t) on the forward path and X{t) = —\9{t)9{T — t) 
on the backward path. Introducing the operator Tr transferring electrons through the 
system in the direction of the current, as well as its counterpart T^, we can write 

Note that = Tl in any system. Consequently, writing out ([7]) explicitly in terms of the 
time-ordered and anti-time-ordered products, one arrives at the conjugation property, 

X*(A)=x(-A). (9) 

We now allow A(t) to be an arbitrary function on the Keldysh contour, X±(t) on the for- 
ward/backward path. Then a generalised expression for the generating function is 

x[A_(t),A+(t)] = (Tce-/c'^*^AW). (10) 
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Assume that the measuring field changes only very slowly in time. Then one can write 



x[A_(t), A+(t)] = exp 



r 



W[A_(t),A+(t)]rft 



(11) 



where W(A_,A+) is the adiahatic potential. Note that Eq. (fTT|) captures the leading, linear 
in T, contribution to the phase. For a truly slow- varying measuring field (analytic in t), 
the corrections are exponentially small in T. Even if the function A(t) has isolated jumps, 
Eq. (|TT|) still holds at large T, but the corrections are then of the order of In T. This is in full 
analogy to the non-equilibrium version 2l|] of the old X-ray problem. Once the adiabatic 
potential is known, the statistics is given by 



lnx(A) = -zrW(A,-A). (12) 
Alternatively one can level off the \± functions in Eq. (fTO!) to different constants as 

x[A_(t),A+(t)]^x(A-,A+), (13) 
then x(A) = x(A, —A) . Note that the conjugation property ([9]) now generalises to 

X*(A„,A+) = x(A+,A_), (14) 



or 



W*(A_,Ah 



-W(A+,A_ 



(15) 



To compute the adiabatic potential we observe that according to the non-equilibrium 
version of the the Feynman-Hellmann theorem 21| , 



d 



■W(A_,A„ 



9A_ 



(16) 



where the averages are defined as 



Tc{A{t)e o 



(17) 



X(A-,A+) 

(and similarly for multi-point functions) where A's are understood to be different on the 
two time branches. Note that the above one-point averages depend on the branch the time 
t is on (though not on the value of t on that branch): 



(A(t_)),^(A(t+)),, 



(18) 
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Therefore the average in Eq. (fT6l) must be taken on the forward branch of the Keldysh 
contour. An advantage of this Hamiltonian approach is that the calculation of the adiabatic 
potential U is reduced to a calculation of Green's functions (GF), albeit non-equilibrium 
ones. So we can use the well developed diagram technique (and relate to many known 
results within this method) without being restricted to the scattering problem. 

The conventional test model to consider is the Anderson impurity model (see 22| for a 
review). The Hamiltonian of the model consists of three contributions, 

H = Ho + Ht + Hc. (19) 

The kinetic part 

Ho = Y, HobpR/L,a] + ^(Ao + ah)dld, , (20) 

cr cr 

describes a single fermionic level (which we shall also call the 'dot') with electron creation 
operators, (ij. (cr is the spin index), energy, Aq, and subject to a local magnetic field, h. The 
electrons in the two non-interacting metallic leads, i = R,L, are represented by ipR/L,a- The 
leads and the dot are coupled via the tunnelling operator, 

Ht = J2 [iLe'^^'^^'di^PLa + iRi^Lda^L + H.c] , (21) 

cr 

with, in general, different amplitudes, 7r,l- We have included the counting field into the 
Hamiltonian (left junction). One can as well incorporate it into the right junction (or both), 
the results are the same due to the gauge symmetry of the Hamiltonian. Finally, we include 
the Coulomb repulsion on the dot. 

He = Un^n^ , (22) 

where n„ = d\dfj. The bias voltage V is incorporated into the full Hamiltonian by assuming 
different chemical potentials in the leads: ^ f^R = > 0. 
It is useful to define two auxiliary GFs, 

Fx{t,t') = -t{Tc{Mt)dKt')})x 

hM = -t{Tc{d{t)ij[{t')})^ . (23) 

Hence the derivative of the adiabatic potential is given by 

d 



X 



W(A_,A+) = ^lim f^e- 

(24) 
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As is standard, the mixed GF's can be written as combinations of bare lead GF's and an 
exact impurity GF, D{t,t'), 



dt"gL{t-t")e-'^^'"^D{t",t') 



FxM = ldt"D{t,t')e-''('"^gL{t"-t'). 



(25) 



Plugging this back into Eg. fl24|] one obtains 



(26) 



where A = A_ — A+. Thus the problem is now reduced to calculation of the impurity GF. 
To illustrate how the method works we start with the trivial case of the resonant level 



model: f/ = 0. Using the GFs of the lead electrons (see, for example, Ref. 16| ) 



9i"{.^) = gt^i.^) = ^27rpoK - 1/2] 
gi'^iuj) = i2TTponi , 



(27) 



where po is the density of states in the electrodes in the vicinity of Fermi level, one easily 
obtains the bare impurity GF (we use the original notation of Keldysh for the GFs because 
the standard identity g + g^^ = g ^ + g^ is spoiled by the measuring field and the 
Keldysh rotation is useless): 



u-Ao- iTL{2nL - 1) - iTR{2nR - 1) 2^6^^/^^^^^ ^ 2iTRnR 

-2^e"*^/2r^(l _ _ 2,r^(l - ur) -uj + Ao- trL{2nL - 1) - trR{2nji - 1) 

where Trl = {.'^Po1r,lY- Inverting this matrix results in 
1 



X 



uj -Ao + iVL{2nL - 1) + iTR{2nR - 1) 226^^/^^^^^ ^ 2iVRnR 

-2ie-'~^I^TL{l - ul) - 2iTr{1 - ur) -c^ + Aq + iTL{2nL - 1) + iTR{2nR - 1) 



where F = F^^ + F^ and we call the object 
V^{uj) = {uj- Aof + F^ + ATlTr U(1 - nB)(e'^/2 - 1) + n^(l - nL){e-'~^'^ - I] 



(28) 
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the 'Keldysh denominator'. On the technical side, while this object is a smooth function 
of energy in the standard technique (expressible via the advanced and retarded GF's), here 
it has discontinuities at the Fermi levels. Inserting this result into Eq. fl26p and integrating 
over A_ (which can be done in a closed form) leads to the formula with the Breit-Wigner 
transmission coefficient, 



as expected for the resonant level problem. 

Turning to the interacting case, we define the impurity self-energy in the standard fashion 

D{u;) = Doicu) + D{iu)±{uj)Do{u;) . (30) 

Substituting this into (l26l) one obtains a general formula for the statistics in interacting 
one-channel systems 

oo 

^ -W(A„, A+) = -r,. [ \2Tr \e'~'/'nL{l - rm) - e-''/'nn{l - ul)] (31) 



-oo 



that expresses the generating function in terms of the (A-dependent) impurity self-energy. 
Here 'Dlu) is the determinant of the (inverse) interacting impurity GF. For A = this 
equation yields the electric current and can be shown to reproduce the result by Meir and 



Wingreen [23|. Formula (13T!) is not restricted to the Anderson model but is applicable for any 
similar one-channel systems (for example, one can add on the electron-phonon interaction 
or consider a double dot). 

The general formula flHT]) allows us to investigate the important limit: linear response 
statistics at zero temperature. Indeed let us take a closer look at the general formula. An 
important technical observation is that the limits V and — do not commute in the 
presence of the counting field. Indeed, calculating the Keldysh determinant in both limits 
we see that even for the non-interacting case 

lim limPo(c^, V,A) = A^ + F^ , (32) 



but 



lim lim Vo{uj, V, A) = + F^ + 4FiFR(e'^ - 1) . (33) 

V—>0 oj— >0 
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The latter scheme needs to be implemented when analyzing the first term in Eg. (1311) in the 
linear response limit, as the energy integration here is restricted to [0, V"]. This leads to 
a transmission coefficient type contribution to the generating function (see below). In the 
second term in Ea. fl3ip. however, the integration over uj is not restricted to [0, V"]. But due 
to Auger type effects 2J] one expects that there are contributions to the current (and PCS) 
at all energies. This effect is itself proportional to the applied voltage and only leads to 
non-linear corrections to the PCS. Hence the energy integration can in fact be regarded as 
restricted to [0, V] even in the second term in Eq. (!3Tl) . Since the self-energy does not have 
external lines and all the internal frequencies have to be integrated over, the limits V ^ Q 
and — s> in this case commute. That means that for the evaluation of the self-energy to 
the lowest order in V one is allowed to use the equilibrium GPs, calculated in presence of 
the counting field A, i. e. f l28|) with ur = ul = np and with the corresponding Keldysh 
denominator. Therefore all diagonal Keldysh GPs are equal to those in the equilibrium and 
all off-diagonal ones are simply proportional to the same diagrams as in equilibrium. Since 
any given off-diagonal self-energy diagram describes an inelastic process, it should vanish 
for a; and we arrive at a conclusion that 

""l 



limE(w) = ReS^(O) 







even at finite A. Eq.f l3ip thus leads to the important result 
Inx(A) = A^ln<^ 1 



(34) 



(35) 



[ReS(^)(0)]2 + r2' 

or to Inx(A) = i\N for the symmetric Anderson impurity model. In the case of an asym- 
metrically coupled impurity, Tr ^ F^, the numerator of fl35l) changes to F/fFi while the 
denominator contains {Tr + Fj^)/2 instead of F. 

The result (1351) allows simple generalisations to asymmetric systems in a magnetic field 
h. According to [251] the real part of the self-energy is given by 

ReEf)(0) =Xc/« + aXs/^, 

where Xc/s are exact charge/spin susceptibilities and fi;~Ao + f//2isa particle-hole sym- 
metry breaking field. Consequently 



Inx(A) 



N 



In 



X 



1 + 



1 + 



F^ 

- (e'> 

[Xc^^-Xshf + T^^ 



(36) 
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An advantage of this formula is that the susceptibihties are known exactly from the Bethe- 



Ansatz results 



26 



271]. We stress that the result fl35|) is not limited to the Anderson model 



but will hold for any similar model, hence the binomial theorem. It is clear in hindsight that 
all the non-elastic processes fall out in the T = linear response limit. Still it is a remarkable 
result that all the moments have a simple expression in terms a single number: the effective 
transmission coefficient. The binomial distribution is universal as long as systems with a 
single conducting channel are concerned. 

III. FCS FOR THE ANDERSON MODEL 

A. Perturbative expansion in the Coulomb interaction 



■ (37) 



We now proceed with the perturbative expansion in the Coulomb interaction U. The 
self-energy, up to f/^-order, is given, in the time domain, by 

~^UD^-{0) + U'[D^~{t)rD^-{-t) -U'[D^^t)]'Dt {-t) 

-U^[Dtit)?D^^-t) tUD^^iO) + U'[D^^it)]'D^^i-t) 

We restrict the calculation to the case of the particle-hole symmetric Anderson model Aq = 
—U/2,. It can be shown that the contribution to the statistics linear in U vanishes in the 
symmetric case. We therefore concentrate now on the second-order correction. 
The equilibrium self-energy is, in fact, known to all orders in U 25l |: 



1 
-1 



2T 



(38) 



sign(ct;) 26{—uj) 
—29{uj) sign(co') 

where Xe/o are the exact even-odd susceptibilities (i.e. correlations of with and n<^ 
respectively), which in weak coupling expand as: 



Xe 



1+3 



TT 



+ 



Xo 



4 J 7r2r2 

For finite V and A in the region —V/2<uj<V/2we find: 



U_ 



(39) 



EM 



^1 - Xc)^ 

iv 



1 
-1 



(40) 



87r2r3 



—e 



-2i\ (ZV 



(f +-) 



QujV 

2 



e--(f -c.)V3(^-c.)^ 
^e-'^{\ + ujY 6ujV 
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Further, for io >V jl one obtains 



-iXjj2 



UJ 



3V 



UJ 



(41) 
(42) 



while for uj < V/2 the following holds: 



[UJ 



UJ 



W_ 



UJ 



(43) 



Svr^rs 

Substituting these self-energies into fl3T|) . expanding around the perfect transmission (hence 
the sign change of A), and formally expressing the result in terms of the susceptibilities, one 
finds the following formula: 



Inx(A) =N {iX + 



BP 



1) + ^(e-2»A _ 



0{V' 



(44) 



This formula is only valid at the order U"^. There are, however, reasons to think that it might 
be exact (see below). It is tempting to interpret the appearance of the double exponential 
terms as an indication of a coherent tunnelling of electron pairs (caution: similar terms would 
also appear for the non-interacting resonant-level model due to the energy dependence of 
the transmission coefficient). 



B. Strong coupling expansion 



In the opposite limit of large [/, the Schrieffer-Wolf type transformation 



tailored to 



the lead geometry 29|], can be applied and results in a Kondo type model. For the latter 
model in the strong-coupling limit, when the spin on the dot is absorbed into the Fermi 
sea forming a singlet, Nozieres 30|] devised a Landau-Fermi-liquid description based on a 
'molecular field' expansion of the phase shift of the s-wave electrons: 



5„{e) = 5q + as + (f)°-{n^ - n^) 



(45) 



where 6q = it/2, a, and 0° are phenomenological parameters corresponding to the residual 
potential scattering and the residual interactions, respectively. These processes are generated 
by polarizing the Kondo singlet and so are of the order ~ where Tk is the Kondo 

temperature. The specific heat coefficient is proportional to a/^iru) while the magnetic 
susceptibility is proportional to the sum a/ {iru) + 2(j)"- /ti. Simple arguments were advanced 
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in Ref. [30j to the effect that, because the Kondo singularity is tied up to the Fermi level, 
there exists a relation a = 2u(j)"' between the two processes in Eq. ( H5|) . In particular, 
this explains why the Wilson ratio is equal to 2. The strong-coupling Hamiltonian that 
describes the scattering and interaction processes encoded in Nozieres Eq. fH^ is of the 
form H = Hq + Hsc + -f^mt- The free Hamiltonian here is 

p,cr 

where is the creation operator for the s-wave electrons, is the creation operator of the 



p-wave electrons, included in order to account for the transport [3l|, and the operator 



p,a 

stands for the (minus) charge transferred across the junction. The scattering term is 

a 

p,p',cr 

while the interaction term reads 



Hi^t = — ^c\c Ac, , (48) 

where = c^ ^ and we have changed to the dimensionless amplitudes a and 0, so that 
in the actual Kondo model a = (p = 1 (in the intermediate calculations it is convenient to 
treat a and as free parameters though). By the nature of the strong-coupling fixed point, 
the operators a and (p are irrelevant in the renormalization group sense and therefore the 
perturbative expansion in a and is expected to converge. 



32 



We now turn to the 



The shot noise in this model was recently discussed in Refs. 
PCS. To this end we introduce the the measuring field, which couples, in the Lagrangian 
formulation, to the current via a term in the action J dt\{t)Q{t) = — J dtX(t)Q(t) that can 
be gauged away by the canonical transformation 

c ^ ca = cos(A/4)c — 2 sin(A/4)a , ^^^^ 
a — > = — z sin(A/4)c + cos(A/4)a . 

We therefore reach the conclusion that the charge measuring field enters this problem as 
a rotation of the strong-coupling basis of the s- and the p-states. While Hq is invariant 
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under this substitution, it should be performed in both the scattering and the interaction 
Hamiltonians, Hsc[c] + -ffint[c] H\ = Hsc[cx] + Hi^t[cx], when calculating the statistics. It 
is easily checked that at the first order in A: = + H^^t + (A/4)/bs + O(A^), where Jbs 
is the backscattering current operator 



p,p ,cr 



- 4Ca)4ca , (50) 



alternatively available from the commutator Jbs = —Q = H]. 

Appl ying the standard linked cluster expansion (still valid on the Keldysh contour, of 
course) 34], we see that the leading correction to the distribution function is given by a 
connected average 

Inx(A) = zNX-l- [ dhdt2{Tc{Hx{ti)Hx{t2)})c + ... (51) 
Jc 

The neglected terms a'^, 0"^, etc., are of the higher order in voltage (temperature) than 

the main correction because of the irrelevant nature of the perturbation. In order to make 
progress with Eq. flSTj) . one only needs the Green's function of the A-rotated c-operator, 
which is easily seen to be the following matrix in Keldysh space: 

gx{p, iv) = zn6{ep - u) { [f{u; - V/2) + f{u; - V/2) - l]fo 
+ [e-'^/^fiu - V/2) + e'^'^f{u + V/2)]f+ 

- [(1 - /(.; - V/2))e^^l^ + (1 - f{uj + r/2))e-^/2]f„ } , (52) 

where fj is the standard choice of Pauli matrices and fi^j) is the Fermi distribution function. 
The correction to the distribution function due to the scattering term P71) is: 

5„lnx(A) = ~ 4^2^2y2 ^i^pi+^P2)^ / dtidt2gp.,{t2,ti)gp^{ti,t2) 

= ^ [ duu^e^'^ - 1)^^(1 - nj,) + (e^^ - l)n^(l - ul)] , (53) 
which, at zero temperature, contributes to Eq. flHT]) a term, 

'5«lnx(A) = y^(e-^^-l). (54) 
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Regarding the correction to the charge distribution coming from the interaction term 
l8|) . any diagrams with a single insertion of the Green's function vanish (therefore there is 
also no acj) cross term) and the only remaining connected graph yields: 

5<^lnx(A) = -^^J-^ j^dhdhg{hMf9it2^h? (55) 



oo 



dt 



K JC 

cos^[\/2 + {yt)/2)] 



TT^Tl J {t + la] 



4 



-oo 



127rT| ^ ' 67rT| 

Combining the results we find that the zero-temperature charge distribution function is 
of the form: 

Inx(A) = ziVA + [{a" + 4>'){e-^^ - 1) 



l2TxTf^ 



,2t^-2iX _^^-]^^Q^y^y (56) 

Let us now try to connect these results to the previous weak coupling calculation. The 
weak coupling expansion of the susceptibilities is given in Section IIII A[ In the strong 
coupling limit we have: 

Xe = {Ta)/{T:TK) Xo = {T<P)/{t:Tk)^ (57) 



where in fact a = (f) = 1 and is the Kondo temperature up to a pre-factor [22|, |26|, |2j 
The programme of extending a Fermi liquid approach to non-equilibrium properties of the 
Anderson model has not been comprehensively carried out yet. There is a Fermi-liquid 



proof, due to Oguri 



that the leading non-equilibrium correction to the zero-temperature 



current is of the form 

4s = Y^^(Xe + 5Xo), (58) 

which is valid for all U and interpolates between the weak-coupling and the strong-coupling 
regimes of the Anderson model. We see that the above result for /bs is simply the strong- 
coupling limit of Oguri's formula. As to the noise and higher moments no analogous Fermi- 
liquid results exist, to the best of our knowledge. However, we would like put forward a 
hypothesis that Eq. (jHj) does represent such a generalisation. Indeed, we see that this formula 
is correct up to the f/^ order in weak coupling, it holds in the strong coupling limit, and it 
reproduces correctly the mean current at all orders in U. 
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FIG. 1: A sketch of the distribution of particles in the leads when > fJ'R, where and fiR are 
the chemical potentials in the two leads. 

IV. USING THE BETHE ANSATZ TO COMPUTE CURRENT NOISE 

In the first part of this paper, we have focused on using Fermi liquid theory to compute the 
leading non-trivial correction (i.e. 0{V^)) of the generating functional, fo^^ both strong 

and weak coupling. We have posited that this computation is exact and have performed 
a number of checks indicating that this is so. In the second part of the paper we adopt 
a different tact, instead focusing upon general features of the moments of the current (in 
particular the noise) in the strong coupling regime over a range of voltages (as measured in 
terms of T^, the Kondo temperature). The tool we use to compute the noise is the Bethe 
ansatz. 

Under the Landauer-Biittiker formalism 



38 



the DC noise, S, reduces to an expression 



given solely in terms of the transmission amplitude, T^-, of electronic excitations across the 
dot: 

^ = tY. d^iT^-Tl). (59) 

Here we imagine the chemical potential of the left lead is greater than that of the right lead 
and so the transmission amplitude, T^-, governs the flow of electronic scattering states from 
the left to the right lead. 

The Bethe ansatz gives us the ability to compute To-(e) in equilibrium up to corrections 
of C(Tx/"\/f/r). We intend to use these equilibrium amplitudes to characterize qualitatively 
the behaviour of the noise as a function of fii — fi^. Using equilibrium scattering ampli- 
tudes necessarily involves missing some of the physics present in a non-equilibrium setting. 
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Nonetheless we believe that this approach yields important insights. In particular, the gen- 
eral behaviour of the noise in a magnetic field that arises from this approach, should be 
robust enough to survive this particular approximation. This is already true of the current 
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with the observed enhancement 
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431] yields results consistent 



through the dot in a magnetic field where this approach 

in the conductance when Hl — IJ^r^ H. 



A. Computation of the Equilibrium Tcr(e) Using the Bethe Ansatz 

We now turn to a brief description of how to compute TL(e) using the Bethe ansatz 
solution of the Anderson model. For additional details see 39[. The description comes in 
three parts. We first describe how to map the equilibrium (zero voltage) problem onto an 
one-channel Anderson model, a model directly solvable by Bethe ansatz. Although we work 
directly in a one lead formulation of the problem, we are still able to make contact with 
scattering in the original two lead picture. We so identify the relevant elements of the exact 
one-channel solution for computing scattering amplitudes. 

To reformulate the problem in a way amenable to the Bethe ansatz, we introduce 
even/odd electrons 

ipe/o = I } „ (tlV^l ± Ir-^r)- (60) 
V 72 + Tr 

With this, the odd electron decouples from the dot leaving us with an interacting theory of 
even electrons alone: 

= ^ y dx^- ii)\^{x)d.^'\\)ea{x) 

+ (7i + ilfi^K^mMd^ + 4v^e.(x)) 

+ Ao^n^ + t/n|n|. (61) 

We point out that the Bethe ansatz could be directly applied to the two channel problem. 
However under the diagonalization of carried out by the Bethe ansatz, the map fl60p is 
implemented implicitly. As such we prefer to make the change of basis directly. 

Under this one-lead reformulation, we are still able to make contact with the scattering 
amplitudes of electronic excitations off the quantum dot. Let T(e)/i?(e) be the transmis- 
sion/refiection amplitudes of electronic excitations of energy, e, between leads in the original 
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two lead picture. On the other hand, the even/odd excitations will scatter off the dot with 
some pure phase, (5e(e)/5o(e). As the scattering in the odd channel is trivial, So{e) = 0. The 
two sets of amplitudes are related according to (pU]l : 



^^Soie) ^ i = R(e)-T{e). (62) 

Our focus will henceforth be on computing Sg. 

To determine (5e(e), we employ an energetics argument of the sort used by N. Andrei in the 



computation of the magnetoresistance of impurities in a bulk metal [40|| . Imagine adding an 
electron to the system. Through periodic boundary conditions, its momentum is quantized, 
p = 27Tn/L. If the dot was absent, the quantization condition would be determined solely by 
the conditions in the bulk of the system and we would write, Pbuik = '27m/ L. Upon including 
the dot, this bulk momentum is shifted by a term scaling as 1/L. The quantization condition 
is then rewritten as 

p = 27™/L = pbuik + 5e(e)/-^^, (63) 

where L is the system's length. The coefficient of the 1/L term is identified with the 
scattering phase of the electron off the dot. We thus must compute the impurity momenta 
of excitations in the problem. 

The Bethe ansatz solution of the one channel Anderson model was first described in 



4l| and [42]. As with any problem with an SU{2) symmetry, the Bethe ansatz yields a 
set of quantization conditions describing two types of excitations, one parameterized by k 
and associated roughly with charge excitations, and one parameterized by A and associated 
approximately with spin excitations: 



11 g{k,) -K-i/ 



2' 



A - g{kj) + i/2 ^ 1^ A^ - A^ + i 

\}^Xc.-9{k,)-i/2 11a.-A^_z' ^ ^ 



(3=1 

where 



5{k) = -2tan-i(r/(fc-erf)); 



g{k) = {k~ea- U/2f/2UT- 
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T = {V^ + V^). (65) 

Here is the total number of particles in the system and M marks out the spin projection 
of the system, 25"^ = iV - 2M (in zero magnetic field M = N/2). 

The possible solutions to the above quantization conditions are manifold. However most 
are only significant at finite temperature. At zero temperature, the ground state of the 
system is formed solely from real fc's (and then only if the magnetic field is non-zero) and 
from bound states of real A's together with complex /c's. Specifically, the ground state 
contains: 

i) N - 2M real k^'s; 

ii) M real A^s; 

iii) and associated with each of the M A^s are 
two complex k's, k^, described by 

g{kl) = gix{Xa) T w{K)) = K± 

x(A) = U/2 + e, - VUT{\ + (A^ + 1/4)^/2)^2. 

y(A) = vW(-A+ (A2 + 1/4)^/2)1/2_ (gg) 

These elementary excitations implement a spin-charge separation in the model with the /c's 
representing the charge sector while the A's represent the spin sector. 

In the continuum limit, these excitations are described by smooth densities, p{k) for the 
real /c's and cr(A) for the A's. Equations valid at the symmetric point of the Anderson model 



describing these densities can be derived in the standard fashion 4l|, |42| : 



P{k) = — + ^ ' 



2n L 

+g'ik) dk'R{g{k)-g{k')){^ + ^ 
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-g'{k) l"" dk'p{k')R{g{k)-g{k'))- 



D 



- [ dkp{k)s{\- g{k)). (67) 

where L is the system size and 

A{k) = dk6{k)/27r, 

s(x) = — -. (68) 

^ ' 2cosh(7rx) ^ ^ 

B marks out the 'Fermi- surface' of the k distribution. Between —D (the bottom of the 
band) and B there is a sea of k excitations. The Fermi surface of the A particles on the 
other hand is set at — oo: at the symmetric point, the sea of A excitations in the ground 
state extends from A = — oo to X = Q, where Q is the bandwidth of the A excitations. This 
is a crucial simplification which makes possible many of our closed-form results. For most 
purposes both bandwidths, D and Q, can be taken to be oo. 

The density equations neatly divide into bulk and impurity pieces via 

p{k) phnikik) + pi^p{k) / L; 

(r{X) ^ fTbuik(A) + c7,^p{X)/L. (69) 

The impurity densities of states contain all the information needed about degrees of freedom 
living on the quantum dot. The equations governing these densities are 



p,^p{k)=A{k)+g\k) / dk'R{g{k) - g{k'))A{k') 

J —oo 

-g'{k) r dk'p,^^{k')R{g{k) - g{k'))- 

J-D 

/oo 
dkA{k)s{X- g{k)) 
-oo 

-B 

dkp,^^{k)s{X-g{k)). (70) 

D 
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For example, the total numbers of spin | and | electrons living on the dot are 



1 1 

2 + 2 



J —oo 

dkp[^p(^k)] 



oo 

B 



oo 



oo 



Udi = / rfA(Timp(A) 



1 1 

2 ~ 2 



B 



dkpimpik). (71) 



The latter equations for each of Udi and Udi are a result of simplifications at the symmetric 
point. 

The energies and momenta of these excitations can be derived through well known tech- 



niques 391]. The energies are given by 



e{k) = k-^-2 I d\ x{\)s{\ - g{k)) 

B 



dk'g'{k')e{k')R{g{k)-g{k'))- 



D 



e(A) = 2x(A) - 2 / d\'R{\ - A')x(A') 



+ / dkg\k)e{k)s{g{k)-X). (72) 

The momenta are akin to the densities in that they divide into bulk and impurity pieces 
391]. The bulk momenta are related directly to the energies via 



e{k) = p{k) - |; 

6(A) = p{X). (73) 
The impurity momenta can be expressed in terms of the impurity density of states 

dkPimp{k) = 2iTpi^p{k); 

dxPimpW = 2Trau^p{\). (74) 
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As already discussed, the impurity momenta are the quantities crucial to computing scat- 
tering phases. These relations will thus allow us to express the scattering phases in terms 
of integrals over the impurity density of states. 

In order to determine the scattering phase of an electron (as opposed to a spin or charge 
excitation), we must specify how to glue together a spin and a charge excitation to form the 
electron. The situation is analogous to adding a single particle excitation in the attractive 
Hubbard model. Adding a single spin f electron to the system demands that we add a 
real k (charge) excitation. But at the same time we create a hole at some A in the spin 
distribution. The number of the available slots in the spin distribution is determined by the 
total number of electrons in the system. Adding an electron to the system thus opens up 
an additional slot in the A-distribution. The electron scattering phase off the impurity is 
then the difference of the right-moving k- impurity momentum, Pimp{k), and the left-moving 
A-hole impurity momentum — Pimp(A): 

= pimp = Pimp(fc) + Pimp(A). (75) 

We must now consider how to choose k and A. 

As we are interested in the DC noise, we must compute scattering away from the Fermi 
surface. Thus if we are to compute the scattering of an excitation of energy, eei, we must 
choose the k and A such that 

6,; = e(A;)-e(A). (76) 

However this constraint does not uniquely specify a particular choice of {k, A). We, in general, 
cannot lift this degeneracy. However at the symmetric point of the Anderson model, we can 



make an ansatz which has already proven to 



De successful in the computation of the finite 



temperature linear response conductance 39|]. The behaviour of the electron scattering 



phase is determined by the impurity densities, pimp and o"imp- At the symmetric point of the 
Anderson model, the scattering phase is expected to vary as ~ T^, the Kondo temperature. 
Of the two impurity densities, only pimp varies as while (Jimp is controlled by the much 
larger scale, VUT. Thus in computing electronic scattering phases away from the Fermi 
surface at T = 0, it is natural to choose A at its Fermi surface value, i.e. A = — oo and vary 
k according to the energy in which we are interested. Specifically, we choose k such that 

eik) = e,i. (77) 
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We thus have removed the ambiguity in the choice of {k, A). 

In making this ansatz we are effectively doing the following. Imagine an electron in the 
leads with some energy Sgi. We can imagine expanding this state in terms of the basis of 
Bethe ansatz states: 

\el>=J2cei,s\s)^ (78) 

s 

where the states \s) are exaci eigenf unctions of the Hamiltonian. Although we only possess 
incomplete knowledge of this expansion, it would be a reasonable guess that it contains 
multiple terms. However our ansatz supposes only a single state contributes. But because 
of the hierarchy of scales, <^ VUT, in the problem, we expect additional terms in the 
expansion of Eqn. flTHl) to have coefficients of 0{Tk/ \/UT). 

Under this ansatz, the scattering phase of the spin ] electron at some energy, e^u above 
the Fermi surface is then 

5l{eei) = Pimpik) +Pimp(A = -oo); 



271 dAaimp(A) + 27r / dk'pimp{k'), 
Jq J-d 



e{k) = e,i. (79) 

When the magnetic field, is 0, e{k) > and we can only directly compute the scattering 
of spin I electrons. However with H > 0, e{k) takes on negative values and so we can also 
compute spin | hole scattering. To add a spin | hole with energy, ehoie > 0, we remove a k 
and a A-hole in the spirit of our previous ansatz: 

e{k) = —ehoie] 

A = -oo. (80) 

The scattering phase is then 

Sloiehole) = Pimp(fc) +Pimp(A) 



2n / (iAa-imp(A) + 2n dk'pi^pik') 
Jq J~d 
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e{k) = -ehoie- 



(81) 



So far we have computed the scattering of spin | objects. To compute spin J, quantities, we 
employ the particle- hole transformation relating spin | to spin |: 



d] — >• df, 



(82) 



where the last line only follows at the symmetric point. With this transformation, we obtain 

S^ei{eei) = Sl^ieho = eel); 



(83) 



Our inabi 



ansatz 



39|. 



ity to directly compute spin | scattering is a technical peculiarity of the Bethe 



B. Nature of the Approximation 

We have now described how to compute the scattering amplitudes as a function of energy 
in equilibrium. The error in using these amplitudes in describing out-of-equilibrium quanti- 
ties such as the current or the noise has two possible sources. The first source can be seen 
from the way the differing chemical potentials of the right and left reservoirs appear in the 
Hamiltonian: 

'Hf, = Hl j dx tpl{x)'4)L{x) + Hr j dx iP]^{x)'4)r{x). (84) 

Under the map to the even and odd basis (!60|) . this term becomes non-diagonal, coupling the 
even and odd sectors. This coupling between sectors, in turn, lifts the model's integrability, 
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the basis on which we compute T'cr(e). One might have hoped 39| that the state of the 
system in non-equihbrium could still be characterized by using in-equilibrium data in a 
fashion analogous to the manner in which out-of-equilibrium quantum Hall edges can be 



characterized exactly 



44l |. While in analyzing the Hall edges, an even-odd transformation of 



bosonic degrees of freedom is employed, the Hall case is much simpler as the current as well 
as the coupling to the voltage bias, can be expressed directly in terms of the odd degree of 
freedom. But in the case of the Anderson model, the current and voltage involve both even 
and odd degrees of freedom. In mapping back to the left-right basis, some sort of breaking 
of integrability then occurs. We will see this explicitly when we compute the noise at zero 
magnetic field and compare it to the exact Fermi liquid results. 

A second source of error in using the equilibrium scattering amplitudes concerns the 
manner in which we construct the electronic scattering states. As we have already indicated, 
there is a multiplicity of choices in how we construct the scattering states. The particular 
choice we employed was the simplest that met the requirement that the energy dependence 
of the scattering varies on the Kondo scale. However this choice is not unique and while 
it produced excellent results for the behaviour of the finite temperature linear response 
conductance, it may not be optimal for the description of out-of-equilibrium scattering. 



V. COMPUTATION OF NOISE IN ZERO FIELD: COMPARISON BETWEEN 
FERMI LIQUID THEORY AND BETHE ANSATZ 

In this section we compute the noise in zero magnetic field and compare it at small 
voltages to the Fermi liquid results from the first part of this article. To compute the noise, 
we imagine biasing the leads as in Figure 1 with /i/j < fi^. For convenience we set fiL to 
zero and fiR = fi. Then as discussed in the previous section, the noise is given by 

S=j iy<T^i^ - ^t) + - Ti)). (85) 

At if = spin I and spin | scattering are the same, i.e. T| = T|. As we are scattering 
electrons from lead R with energy e < into lead L, we are equivalently interested in 
computing the scattering phase of a hole. In Section IV we demonstrated that at = we 
are able to compute the scattering amplitude of a spin | hole. Specifically we have 
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j-Q pk 

Sl^ie) = 2tx / rfAaimp(A) + 2ti / dk' p;^^{k') 

jQ J-D 



e{k) = e. (86) 



In 391 we were able to evaluate these expressions in closed form: 



5T (,) ^ 3^_sin-i(L^!/S) 

OO 

^ 1 + 2n^V2f7r^ 



X 



^^e-^9(fc)(l+2n)j^gj^^-^^|. 



TT TT V 2 

The last equation gives the crossover scale, T^, the Kondo temperature, in terms of the bare 
parameters of the model [45|. With this, Ty^ equals 



Tyi = ^-^ + 0{^^). (88) 



For typical realization of Kondo physics in quantum dots, i.e. 36|, |37|], the error term is 
insignificant. And so we compute the noise to be 



o2 



5(/.) = \ - f, tan-^(l-) ) . (89) 



The quantity S/V is solely a function of the ratio V/Tk and is plotted against in Fig. 

2. (Note that = eV.) We see that as V/T^ is varied, the noise rapidly rises at first, peaks 
at eV ~ — l.lSTfc, and then begins to gradually decline. 

This behaviour is closely related to the Kondo resonance in the spectral impurity density 
of states. As we express the scattering phases in terms of the impurity density of states, we 
probe the resonance. The width, w, of this resonance, w ~ Tk, corresponds to the energy 
scale at which we expect maximal noise. 

The noise in this case can be reexpressed in terms of the current, /, and the differential 
conductance, G, 

J(/i) = -2fr,tan-i(|-); 

J-k 
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with the result 



At small fi, 



S{fi) = + (91) 



Sif^ < 0) = -^^f^'Tl 



2e2 87r2 

"/T'ge 

We can compare this result with the Fermi liquid result 



In making this comparison, there is a certain arbitrariness in how one defines Tk- This can 
be overcome by appealing to the finite temperature linear response conductance, G(T) to fix 
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the manner in which is to be defined. In our conventions then, G{T) = 2e^ /h{l — 
We see then that the Fermi liquid result differs from the result based upon the equilibrium 
Bethe ansatz scattering states by a factor of 5/8. 

A similar difference can be found between the Fermi liquid and Bethe ansatz computation 
for the current. At small /x, the leading order correction to the current (of 0{ix^)) is given 
by 

(92) 

This compares to the Fermi liquid result 

- 'tM- 

We see that the Fermi liquid result is considerably smaller than that of the Bethe ansatz. 
This might well reflect the role of incoherent scattering processes that would be unaccounted 
for properly by using equilibrium scattering amplitudes. 

Finally we consider the value of the effective charge, e*, in the problem. This charge is 
given as a ratio of the noise to the backscattering current: 

e* = S/hs. (94) 



In the case of the Bethe ansatz, we flnd e* = e, that is, we flnd the Johnson-Nyquist result 
for shot noise in the weak scattering limit. However in Fermi liquid theory, the effective 
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FIG. 2: Plot describing the evolution of both the noise, S/V , and the current, I /V , as a function of 
the applied voltage. The scale on the l.h.s. governs the noise while the scale on the r.h.s. governs 
the current. 

charge is found to be e*^^ = 5/3e 321]. That the effective charge goes unchanged from its 
non-interacting value is again presumably a consequence of the use of equilibrium scattering 
states. 

While the use of the equilibrium Bethe ansatz scattering states gets quantitative details 
of the noise (and the current) incorrect, it gets qualitative features correct. In particular the 
conductance as a function of voltage is Lorenztian-like with a width of order Tk- Similarly 
the zero field noise as a function of voltage increases rapidly (on a scale of order Tk) and 
thereafter decreases slowly. Both of these features would be expected to be present based 
solely upon the presence of the Kondo/Abrikosov-Suhl resonance whose width is governed 
by the scale Tk- 



VI. COMPUTATION OF NOISE IN FINITE FIELD 



In this section we compute the noise in a finite magnetic (Zeeman) field. We will here 
argue that the finite field noise possesses features that can be considered a "smoking gun" 
461] of Kondo physics. We believe that these features are robust and so should be captured 
by our approach. 

To compute the noise, 

S=^^f^ de{T^{l - T^) + Ti(l - TJ), (95) 
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in a finite field, we must consider the contributions of the spin | and spin | currents indi- 
vidually. From Section III, T| and T| are given by 



I = sm 



rQ pk 
2ti I rfAo-imp(A) + 27r / dk'pi^p{k') 

JQ J-D 

e{k) = -e; 



n = sin^(^); 



rQ rk 
271 dAcrimp(A) + 27r / dk'pimp{k'), 
Jq J~d 

e{k) = e. (96) 



In 



391 we evaluated these expressions in two cases, H '^Tk and H > Tk 



case i: H <^Tk 

We found for spin | hole scattering 



Slie,^ > 0) = -TT - sin-i i ^-i^'^o-Hf/Tj . 



--sin-i(i— (97) 



while for spin | hole scattering we arrived at 



. . 5 . _i/l-(e,„ + /7/2)Vr2 



4 Vl + (e.o + if/2)VT| 



1 _i fl-HViAff 
--sin ' 



2 Vl + i/7(4T,2) 
Consequently the transmission amplitudes in this case equal 



(98) 
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2 I (1 + m/AflYI\l + {^^, - H/2f/fl) 



(99) 



case ii: H > Tu 



Using a Weiner-Hopf analysis, and were determined in Ref. 39 1 to be 

5l = 7r + 2tan-i(2(/-i-^?(A;))); 



37r 

5l = y +tan-i(2(/-i-6)), 



(100) 



where / ^ sets the Kondo scale, Tk ^ e '^^ \ and / ^ — 6 is given in terms of the ratio H/Tk- 

r-i u r . 



2f/' 



/-^-& = -iog(— ,/^). 



(101) 



/c is parameterized in terms of the energy, e, by the expression 



e{k) = -H(1- :^tan~i 1 



27r g{,k)-h 

1 1 
47r2 1 + {g{k) - by 



;^ + l -(.(.)-. )tan-(-^) 



-C + -log{A7r\l + {gik)-bn 



VWU f 1 e-^^ 



+ e"""^'^^ tan 



^^(fc)-6' 



+0((<?(A;)-6)"3) 
where C = .577216 ... is Euler's constant and b is given by 



b = — log( — \/ — — 
TT H \ ne ' 



With this 



l + (2(/-l-^(A;))2) 



+ ■ 



(102) 



(103) 
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FIG. 3: Plots describing the behaviour of the noise, S/V ^ as a function of the apphed voltage, V ^ 
for a variety of magnetic field values. 

In writing fll04p we have omitted writing terms arising from the full expression for e{k) in 
(11021) . But because of the logarithmic dependence upon H/Tk, such terms are needed if we 
are to compute the noise with reasonable accuracy for fields, H, not far in excess of T^. 
Notice that the spin | scattering does not vary as a function of energy, an approximation 
valid for H ^ Tk. 

We plot the noise, 5/1^1, in Figure 3 for a variety of values of H/Tk- For small H/Tk 
the noise is smooth and without structure. However as we vary H from H < to H > T^, 
a marked variation occurs as seen in top panel of Figure 3. The noise rapidly increases, 
achieving its maximal value of roughly i/ ~ T^, the crossover scale, before again decreasing. 
This behaviour is repeated in the differential noise —dyS, plotted in Figure 4. The noise's 
maximal value at = is a reflection of the maximum in T^(l — T^). Such a maximum 
occurs for T'^ ~ 1/2. Thus at H = Tk, the average transmission amplitude for spin f 
excitations has been reduced roughly by 1/2. 
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At large H/Tk, the noise develops a double humped structure near e\V\ ~ H. This feature 
is more apparent when we examine the differential noise in the lower panel of Figure 4. In 
varying V near the peak, the differential conductance G oc T'^{V), passes through the value 
Gmaxf^ = e^/{2h) twice. As such the quantity, dyS ~ T\l — T^), possesses two peaks. 

Given the bias at which it occurs, the doubled peak is intimately related to the peak in 
the differential conductance seen near eV ~ H. The peak in the differential conductance 



owes its origin to a field induced bifurcation in the Kondo resonance 



46|: as the Kondo 



resonance shifts so does the peak in the conductance. In Refs. [39|, |43|] this bifurcation was 
studied where it was found that the peak occurs at a value of eV distinctly smaller than 
H and not eV = H. In the case of the noise, we again find that the peaked structure in it 
occurs at values of eV smaller than H. 

We believe this double peaked structure in the noise, inasmuch as it depends on the gross 
dependency of the scattering amplitudes upon H, to be a robust feature. Less certain are the 
quantitative predictions that arise from this analysis. Nonetheless we will proceed to analyze 
the structure of the differential noise peaks. We do point out that we have some confidence 
that this analysis has merit as a subset of its corresponding predictions for the behavior of 



47|. In 



the current in a magnetic field have been shown to be at least qualitatively correct 
particular, the prediction that the peak in the differential conductance occurs for values of 
eV smaller than H has been observed in experiments on carbon nanotubes \^^ - 
The differential noise, —dyS, is given by 

dyS = j{t\1-T^)+T^{1-T^)). (105) 

For V ^ H and H ^ T^:, ^ T|, the locations of the peaks again occur when 

= 1/2. 

As = sin^(5^^/2), the scattering phases that correspond to this amplitude are, = 
|/^. From fllOOl) . this in turn implies g{k) = I^^ ±1/2. Using [TIM the biases, V±, at which 
the two peaks occur equal 

H 1 

eV± = (106) 

The height of the peaks can also be determined. The maxima of the peaks occur when 
T^^ = 1/2. Consequently, the height of the peaks in —dyS are given by adding e^/{4h), the 
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FIG. 4: Plots describing the behaviour of the differential noise, -dyS, as a function of the applied 
voltage, V, for a variety of magnetic field values. 

contribution due to the spin | current, to the contribution from spin | scattering with the 
result, 

"^^^ -AhA{i-^-by + i 

This result holds for either peak, a consequence of the lack of variation in for eV, H ^ T^. 

We are also able to compute the full width at half maximum (FWHM) of the peaks. As 
the peak maxima occur for T\l — T^) = 1/4, or phases, Sf^l"'''^ = f ^^e FWHM occur 
for T\1 — T^) = 1/8. For the peak, 5^^"'^^ = 7r/2, the phases of the FWHM then correspond 
to 

ho ~ 4' 4 ■ 

These two phases occur at energies parameterized by values of k given by 
Hence the width of this peak is 

H ( 

e/W = — I tan"^ 



27rV -6+ |tan(|)^ 

-tan-V— \ 7^))- (109) 

V-i-6 + itan(f )^ ' ^ ' 
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FIG. 5: Plots describing the evolution of the differential noise peaks with increasing magnetic field. 
In the top panel are plots of the widths of the two peaks, in the middle panel plots of the two 
peaks' locations, and in the bottom panel, a single plot of both peaks' (identical) height. 

Similarly, for the peak corresponding to S^^'^^ = the phases of the half-maxima are 

cFWHMi _ Svr Ttt 
"ho ~ 4 ' 4 ' 

and consequently, the FWHM of this peak is 

H / 1 
eAV = — I tan"^^ ^ 



27rV V-i -6- itan(|^)' 
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The width of the two peaks is thus notably different, with the peak occurring nearer to 
e\V\ = H the narrower. 

The various peak parameters are plotted in Figure 5. We see in the bottom panel that 
the height of the peaks approaches an asymptotic value of /Ah logarithmically in H/Tk- 
This limit corresponds to a situation where only spin | electrons contribute to the current. 
In the middle panel of Figure 5 are plotted the biases, V±, at which the two peaks occur. In 
the large H/Tk limit, eV± approaches H. However even at large H/Tk there is a significant 
deviation from H, a consequence of the logarithmic dependence upon H/T^. This behaviour 
mimics that of the peak in the differential conductance: for large H/Tk, this peak occurs at 
a value of eV distinctly smaller than H 39|]. Finally in the top panel of Figure 5 is plotted 
the width of the two peaks. Again the peak that occurs at a bias closer to -H is the more 
narrow of the two. 



VII. CONCLUSIONS 



To summarize, we have presented analyzes of the moments of the current using both 
Fermi liquid perturbation theory and the Bethe ansatz. Both approaches however remain 
incomplete. In the Fermi liquid approach, while we have derived results for the higher 
moments of the current valid, we believe, at all orders of the interaction strength, we have 
done so only at the lowest non-trivial order in the voltage. To expand these computations 
to higher orders in the voltage remains a challenging problem for future research. In terms 
of the Bethe ansatz treatment, we have managed to identify various robust features of the 
current noise that should be experimentally identifiable in current realizations of quantum 
dots. In particular, we have identified a double peaked structure in the noise that appears 
at finite magnetic fields. But because of our use of equilibrium scattering amplitudes, we 
have been able to focus only upon qualitative aspects of the physics. Quantitatively, the 
methodology produces results at variance with Fermi liquid theory. It thus remains an 
important goal to pinpoint the precise origin of this disagreement. 
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